The Kondo effect in the presence of the Rashba spin-orbit interaction 
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We study the temperature scale of the Kondo screening of a magnetic impurity which hybridizes with a two- 
dimensional electron gas in the presence of the Rashba spin-orbit interaction. The problem is mapped to an 
effective single-band impurity model with a hybridization function having an inverse-square-root divergence at 
the bottom of the band. We study the effect of this divergence on the Kondo screening. The problem is solved 
numerically without further approximations using the numerical renormalization group technique. We find that 
the Rashba interaction leads to a small variation of the Kondo temperature (increase or decrease) which depends 
on the values of the impurity parameters. 
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The spin-orbit (SO) interaction is a relativistic effect due 
to the inter-dependence between electric and magnetic fields 
when considered from different reference frames. It leads to 
a coupling between electron's spin and its (orbital) motion in 
real space'. The effect is stronger in heavy elements from 
the bottom of the periodic system. The spin-orbit interaction 
plays a central role in many proposals for spintronic devices^ 
and it is the physical origin of the recently discovered topo- 
logically non-trivial insulator phases which do not break the 
time-reversal symmetry^^. The spin-orbit interaction does 
not break the Kramers degeneracy, thus the Kondo screening 
of magnetic moments is possible^r— and should be observable 
in magnetic adatoms on topological insulator surfaces or on 
thin layers of heavy elements^. 

The Kondo effect due to magnetic impurities in bulk sim- 
ple metals was shown to be strongly suppressed by doping 
the host material with small ammounts of Pt impurities which 
have a large value of the SO coupling constant", although 
these results have been questioned and the expected trend 
with the increasing strength of the the SO coupling has not 
been confirmed in later experiments'". An Anderson impu- 
rity model with SO scattering term for the conduction band 
electrons was studied by performing a Schrieffer- Wolff trans- 
formation into an effective Kondo modeU^ and it captured 
some (but not all) of the observations of Ref. Kondo 
screening has also been studied through weak-localization ef- 
fects; the results indicate that adding spin-orbit scatterers does 
not change the magnetic scattering at all^ii. It was pointed 
out that the spin-orbit scatterers play the same role as elas- 
tic nonmagnetic impurities because they do not break the 
time-reversal symmetry^,, thus the Kondo temperature is ex- 
pected to remain unchanged, in agreement with the exper- 
iment. For similar reason of time-reversal invariance, the 
Kondo screening of magnetic moments is expected to occur 
for magnetic impurities adsorbed on the surfaces of three- 
dimensional topological insulators-^, unless the doping level 
is sufficiently high for the spontaneous breaking of the time- 
reversal symmetry due to magnetic ordering. The effect of 
the spin-orbit coupling in a two-dimensional (2D) electron 
gas with Rashba interaction' on the Kondo screening has been 
studied in Ref.[l5l The work was based on a Kondo impurity 
Hamiltonian and after considering to which conduction-band 



angular modes the impurity spin couples to, a two-channel 
Kondo model was derived, in which each channel has a dif- 
ferent dispersion relation. By expanding the dispersion re- 
lations to linear order, remarking that the Fermi velocity is 
the same in both helicity channels, and after some further 
manipulations, the model maps to a single-channel single- 
impurity Kondo model with unchanged dimensionless Kondo 
coupling constant pj. In other words, it was predicted that 
turning on the Rashba coupling (all other parameters remain- 
ing the same) does not change the Kondo temperature at all, 
except perhaps for some effective band-width effects. The 
same problem has been considered starting from an Ander- 
son impurity model and using the Varma-Yafet variational 
ansatz'.^; this work did not specifically address the question 
of the Kondo temperature beyond noting that it remains con- 
stant in the limit of small Rashba coupling, but focused instead 
on the spatial distribution of the spin correlations. Recently, 
a new study of this system used a Schrieffer- Wolff transfor- 
mation to obtain an effective two-channel Kondo model with 
an additional Dzyaloshinsky-Moriya (DM) interaction' Us- 
ing a scaling renormalization group analysis, it was argued 
that the DM term renormalizes the Kondo coupling and pro- 
duces an exponential enhancement of the Kondo temperature. 
This effect has not been observed in previous studies. Since 
all the cited works are based on various mappings and ap- 
proximations, it is unclear which result for the variation of the 
Kondo temperature (or lack thereof) is correct and what ap- 
proximations are responsible for the disagreement. For this 
reason, we reexamine the problem using an approach where 
all the approximations are well controlled. To wit, we inte- 
grate out the conduction-band degrees of freedom to obtain 
an impurity action with a hybridization function which has an 
inverse-square-root-divergence at the bottom of the conduc- 
tion band. This step is exact. We then numerically study the 
effective impurity model using the numerical renormalization 
group (NRG) method'^ii^. While the NRG approach is not ex- 
act, the approximations involved are all well controlled. We 
find that the Kondo temperature is only weakly affected by the 
Rashba interaction: it may weakly (approximately linearly) 
increase or decrease depending on the values of the impurity 
parameters. 

We study a magnetic impurity embedded in the 2D electron 



2 



gas with Rashba spin-orbit interactio n'''^~^^'^° . The impurity 
is described using the single-orbital Anderson impurity model 
Hq with an additional Rashba interaction term i?so" 



Ho = €{n^ + ni) + Uni^ni 



Hso = a^i/'j^ (fcjjCTj, - feyCTj,) V'k 
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The operator rf^ creates an electron in the impurity level, 
while cj^^ correspond to the conduction-band electrons with 
the dispersion = k"^ /2m* + Eq, where k is the crystal mo- 
mentum, m* the effective electron mass, and the bottom of 
the conduction band. The occupancy operator is ria ~ d^d^. 
The chemical potential is set to the energy zero, /i = 0. In 
Hso, parametrizes the strength of the Rashba spin-orbit in- 
teraction, xjjy is a spinor field {ckt, Ckj.}"'", while (Tx/y are 
Pauli matrices, k and (/)k are the polar coordinates of the wave 
number k, with the polar axis oriented along ky. 

We now switch to a continuum representation in a box of 
volume V = 1. The sums are transformed in the standard 
way as — > V/(27r)^ J fcdfcd(/). The operators are trans- 
formed as Ckcr — > (27r/V)'0kcr; they are normalized such that 
Wkcr.^Iv'} - _ k')(5,,.. Noting that 5(2) (k) = 

5{k)5{(j)) I k, we switch to a polar representation using -^kcr — >■ 

l/Vk^k<t,a, with Ufc0<.,Cl'0V'} ^(^' - - (p')Saa'- 

Finally, we switch to the angular momentum basis ^^kipa = 
e*™'^c)!^, where m is the orbital magnetic quantum 

numbe r^' '^''^ and the anticommutation relations take the form 
{Cfc^>(c"l')^} = S{k - k')6„,r,r'Saa'- Collectively, these 
transformation steps can be written as 
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We then introduce the linear combination which diagonalizes 
the conduction-band Hamiltonian including the Rashba cou- 
pling term a'^i'^ 
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where h = ±1 is the chirality quantum number. The corre- 
sponding band energies are 



e/c/i = fife + akh = 



(fc + hkp)^ 
2m* 



Eq — Er, 



(7) 



where fco = m*a is the Rashba momentum, and Er = 
k'^/2m* is the Rashba energy. This implies that the bottom of 
the conduction band shifts from Eq to E^ — Er upon switch- 
ing on the Rashba interaction. 



Figure 1 : (Color online) The dispersion and the density of states in 
the presence of the Rashba spin-orbit coupling. The label ft = ±1 
indicates the helicity of the band. 



We obtain the effective impurity model 
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We now integrate out the conduction-band modes and obtain 
an effective impurity action of the following form: 
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The hybridization function A/j (r) is the Fourier transform of 



1 Z""" 

— / fcdfc— 
47r Jo zw, 



f-kh 



(10) 



While A+(icl'„) and A_(iw„) are somewhat complicated 
functions, their sum A = A^ is simple. Assuming 
Vk = V", the imaginary part r(i;) = -l/(27r)ImA(£; + iO+) 
is found to be of the same functional form as the density of 
states (DOS) of the conduction band: 



E<Eo-Er 

r(^) = ro<| ^ E-(E:-En) Eq-Er<E<Eo (11) 

1 E> Eo 



where To = ^i^^^- This result is remarkably simple and 
could have been guessed in advance, since the impurity is as- 
sumed to be point-like and simply couples to the local DOS 
at its position. The only effect of the Rashba interaction as 
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far as the local impurity properties are concerned is the emer- 
gence of an additional energy interval [Eq — En; Eq] with 
finite density of states which diverges with an inverse-square- 
root divergence at the lower limit, see Fig.[T] The effect of this 
diverging density of states is difficult to evaluate analytically, 
therefore we resort to using a numerical technique. 

We compute the impurity spectral function using the nu- 
merical renormalization group (NRG) methodi^ii^ with ex- 
tensions for arbitrary density of states^^. NRG can han- 
dle diverging DOS both at the Fermi level?^'^-, as well as 
away from the Fermi leveP". Here the discretization has been 
performed using the discretization scheme from Refs.l27ll2i 
this scheme easily handles inverse square root divergencies 
at finite frequencies without a need for introducing artificial 
cut-offs. NRG parameters were A = 2, twist-averaging over 
= 64 discretization grids, and the truncation cutoff set at 
IOwat, where cjjy is the characteristic energy scale at the A^-th 
NRG step. 
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Figure 2: (Color online) The impurity spectral functions for a range 
of the Rashba spin-orbit coupling strengths. The gray curve (red 
online) in the left-hand panels is the effective hybridization function 
that the impurity state mixes with. The parameters are U — 1, e — 
—0.5, and the hybridization strength at Fermi level for En = has 
been set to To = 0.05. The band bottom is Eo = —0.2. 

The calculated spectral functions are shown in Fig.|2l Left 
panels show the spectral function over the full energy interval, 
while the right panels are a close-up on the Kondo resonance 
at the Fermi level. The hybridization function entering the 
impurity model is shown as a gray curve (red online) in the 
left-hand panels. The spectral function features two atomic 
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Figure 3: (Color online) The Kondo temperature calculated using 
the NRG. The parameters are specified in the plot. The temperatures 
are rescaled by the Kondo temperature in the absence of the Rashba 
interaction. 



resonances, one at w = e (corresponding to the extraction of 
an electron) and one atoj^e + U (corresponding to the ad- 
dition of an electron). The width of atomic spectral peaks is 
approximately 2r(a;), where ui is the peak positio n-^ The 
peak at w = e lies outside the band for all values of En con- 
sidered, thus it should have zero width (i.e., it is a delta peak); 
the finite width in the presented spectral functions is due to 
the spectral broadening in the NRG. In addition to the atomic 
peaks, the spectral function features the many-body Kondo 
resonance which peaks near the Fermi level. The long loga- 
rithmic tails of the Kondo resonance are asymmetric due to the 
particular energy dependence of the hybridization function. 
At the bottom of the band, at lj ~ Eq — E^, the spectral func- 
tion has a maximum, then it drops to zero at lower frequencies 
since in that range there are no states in the band. Again, the 
drop to zero in presented spectral functions is overbroadened 
due to technical reasons. 

The Kondo temperature (defined as in Ref. [33l) has been 
determined from the impurity magnetic susceptibihty, but it 
could alternatively be extracted from the Kondo resonance 
peak widths. The results for a several different choices of 
the impurity parameters are shown in Fig. [3] Note that we 
have included parameter sets both with e + U/2 = and with 
e + U/2 0. We find that the Kondo temperature exhibits 
some variation as a function of E^, however the variation is 
not exponential'^, but rather linear in the laige-E^ limit with 
more complex variation for small values of Ej^ (which can be 
attributed to the emergence of the inverse-square-root diver- 
gence in the hybridization function). Depending on the value 
of e, U, and Fq, the Kondo temperature is either an increasing 
or a decreasing function of En- 

In systems with constant hybridization function F, the 
Kondo temperature is given approximately by 



Tk ^ Dy^pJ cxp ( - 
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where pj = SF/ttJ/ is the dimensionless Kondo coupling 
constant. If the hybridization function is energy-dependent, 
we have to use its value at the Fermi level (Fq) in the expres- 
sion above. In the presence of the Rashba interaction, if the 
bottom of the conduction band is initially below the Fermi 
level (i.e., i?o < 0), the density of states at cj = will re- 
main unchanged when the Rashba coupling is increased. This 
explains the absence of an exponential renormalization of Tk 
and is in agreement with Ref. [isl Nevertheless, in Eq. ( fT2] i 
the effective bandwidth D does depend on the details of the 
hybridization function, thus some variation of is in fact 
expected and indeed observed, as shown in Fig. [3] 

We note that the effective impurity action, Eq. (|9]l, has the 
form of a single-band problem. In fact, the single-orbital 
Anderson impurity model is always effectively a single-band 
problem, no matter what kind of conduction-band it couples 
jq8j24 rpjjg required unique effective conduction band can 
always be constructed by taking the combination of states 
which couples to the impurity orbital (here a state propor- 
tional to the linear combination ^kc\^„) as an initial state 
in the Gram-Schmidt orthogonalization procedure applied on 
the conduction-band Hamiltonian. Consequently, one ob- 
tains a single-band tight-binding Hamiltonian, while all other 
conduction-band states are fully decoupled from the impu- 
rity. This argument holds universally. A multi-channel Kondo 
model can be in some cases obtained from the single-orbital 
Anderson impurity model by using a non-optimal basis for 
the conduction-band modes that also includes states which 
are in reality fully decoupled from the impurity. These states 



should not in any way affect impurity properties. Unfortu- 
nately, the complete irrelevance of the redundant conduction- 
band states can only be observed if the problem is solved ex- 
actly. If the problem is, however, approached with an approx- 
imate method, there exist a potential pitfall: the presence of 
the redundant states can affect the impurity properties in some 
spurious way. 

We also note that the problem of an impurity in the two- 
dimensional electron gas can never be truly particle-hole sym- 
metric, because the conduction-band itself is not particle-hole 
symmetric. We therefore do not expect any particular differ- 
ence in the role of the Rashba spin-orbit interaction depending 
on whether e + J7/2 is zero or not'-. 

The predictions of this work could be experimentally tested 
in a system where the spin-orbit interaction can be tuned, 
such as 2D electron systems in semiconductors with metal 
gates. Furthermore, the rapid advances in the field of ultra- 
cold atom system suggest that it might become possible to 
build 2D fermionic gases with tunable SO interaction^— and 
couple it to a magnetic impurity^. 
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